Exact solutions of a Flat Full Causal Bulk viscous FRW cosmological model 

through factorization 



O. Cornejo-Perez 
Facultad de Ingenieria, Universidad Autonoma de Queretaro, 
Centra Universitario Cerro de las Campanas, 76010 Santiago de Queretaro, Mexico 

J. A. Belinchon 

Departamento de Fisica Atomica, Molecular y Nuclear. Universidad Complutense de Madrid, E-28040 Madrid, Espana 

(Date text: January 25, 2013) 

We study the classical flat full causal bulk viscous FRW cosmological model through the factor- 
ization method. The method shows that there exists a relationship between the viscosity parameter 
s and the parameter 7 entering the equations of state of the model. Also, the factorization method 
. allows to find some new exact parametric solutions for difi'erent values of the viscous parameter s. 

y—^ ' Special attention is given to the well known case s = 1/2, for which the cosmological model admits 

I scaling symmetries. Furthermore, some exact parametric solutions for s = 1/2 are obtained through 

. the Lie group method. 
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^ ■ I. INTRODUCTION. 

(N : 

' Factorization of linear second order differential equations is a well established method to find exact solutions through 

^ , algebraic procedures. It was widely used in quantum mechanics and developed since Schrodinger's works on the 
Q^' factorization of the Sturm-Liouville equation. At the present time, very good informative reviews on the factorization 
method can be found in open literature (see for instance [U, 01). However, in recent times the factorization method 
{^JT) has been applied to find exact solutions of nonlinear ordinary differential equations (ODE) 0-0|- based on 

|, previous Berkovich's works Q , it has been provided a systematic way to apply the factorization method to nonlinear 
• second order ODE. In [5|, Wang and Li extended the application to more complex nonlinear second and third order 
^ ^ ODE. The factorization of some ODE may be restricted due to constraints which appear in a natural way within 
QQ . the factorization procedure. However, here it is shown that by performing transformation of coordinates, one can 
' be able to get exact parametric solutions of an ODE which does not allow its factorization or presents cumbersome 
0^ , constraints. 

' The purpose of the present work is to apply the factorization method to study the full causal bulk viscous cosmolog- 
\ ical model with flat FRW symmetries. Since the Misner j8|] suggestion stressing the fact that the observed large scale 
' isotropy of the Universe may be due to the action of the neutrino viscosity when the Universe was about one second 
04 ■ old, there have been numerous works pointing out the importance of the physical processes involving viscous effects 
> in the evolution of the Universe (see for instance [91] ) . Due to such assumption, dissipative processes are supposed to 
play a fundamental role in the evolution of the early Universe. 

The theory of relativistic dissipative fluids, created by Eckart [Tol| and Landau and Lifshitz [Tl'] has many drawbacks, 
and it is known that it is incorrect in several respects mainly those concerning causality and stability. Israel [T2j 
\ formulates a new theory in order to solve these drawbacks. This theory was latter developed by Israel and Stewart 
into what is called transient or extended irreversible thermodynamics. The best currently available theory for 
analyzing dissipative processes in the Universe is the full causal thermodynamics developed by Israel and Stewart 
[T^ . Hiscock and Lindblom [31 and Hiscock and Salmonson [l5|. The full causal bulk viscous thermodynamics has 
been extensively used to study the evolution of the early Universe and some astrophysical process [H, [ol . 

The paper is organized as follows. In Section II, we start by reviewing the main components of a flat bulk viscous 
FRW cosmological model, and introduce the factorization technique as applied to the cosmological model. Field 
equations (FE) of the classical bulk viscous FRW cosmological model [17] reduce to a single nonlinear second order 
ODE, the fundamental dynamical equation for the Hubble rate. By performing a transformation of both the dependent 
and independent variables and using the factorization method, this equation is transformed into a nonlinear first order 
ODE. The order reduction of the equation for the Hubble rate allows to find a variety of new exact parametric solutions 
of the FE for the viscous FRW cosmological model. Furthermore, the factorization technique provides relationships 
for parameters entering the factorized equation. Then, a noteworthy result is that the viscosity parameter s is not 
longer assumed to be independent of the values of parameter 7. Such parameter relationships have not been previously 
reported. In Section HI, several particular models for s ^ 1/2 are studied. We obtain new exact parametric solutions 
through factorization and compare with the ones obtained by several authors [T8l - [25| who use different approaches. 
Section IV is devoted to the special case s — 1/2, for which the model admits scaling symmetries. The scaling solution. 
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previously studied by many authors is obtained. In order to obtain more new solutions and compare the solutions 
obtained through factorization for s — 1/2, we consider the Lie group method for this special case in Section V. Some 
conclusions end up the paper in Section VI. 



II. THE MODEL. 



We consider a flat FRW Universe with line element 

ds^ = -dt^ + f{t){dx^ + dy^ + dz^), (1) 

where the energy-momentum tensor of a bulk viscous cosmological fluid is given by [l7j : 

Tt ^{p + p + Il)u,u'' + {p + n)dt (2) 

where p is the energy density, p the thermodynamic pressure, 11 the bulk viscous pressure and Ui the four-velocity 
satisfying the condition mu^ — —1. We use the units 8ttG = c — 1. The gravitational field equations together with 
the continuity equation, T^^j, — 0, are given as follows 

2H + 3i^2 = - n, (3) 

3i/2 ^ (4) 
Tl + Ttl = -SCH-^Tu(^3H+^-i-^y (5) 

p = -3{^p + n)H, (6) 

where H — f / f. In order to close the system of equations we are assuming the following equations of state p7j 

p=(7-l)p, C = ^ = T^(p-'^ap^-\ (7) 

where T is the temperature, f the bulk viscosity coefficient and r the relaxation time. The parameters satisfy 76 [1,2], 
s > 0, and r = ^1 — . The growth of entropy has the following behavior 

S (t) « ~ikg^ [ UHfT-^dt. (8) 

Jto 

The Israel-Stewart-Hiscock theory is derived under the assumption that the thermodynamical state of the fluid 
is close to equilibrium, i.e., the non-equilibrium bulk viscous pressure should be small when compared to the local 
equilibrium pressure |n| << p = {7 — l)p. Then, we may define the l(t) parameter as: / — |n|/p. If this condition is 
violated then one is effectively assuming that the linear theory also holds in the nonlinear regime far from equilibrium. 
For a fluid description of the matter, the condition ought to be satisfled. 

To see if a cosmological model inflates or not it is convenient to introduce the deceleration parameter q = dH^^ / dt — 
1. The positive sign of the deceleration parameter corresponds to standard decelerating models, whereas the negative 
sign indicates inflation. 

The fundamental dynamical equation for the Hubble rate is given by [T^ 

H - A— + {3H + CH^-^') H + DH^ + EH'^-^" = 0, (9) 

where 

A = {l + r)=2--, B = 3, C^?}-", D = -h-2), E^-3^-'-f. (10) 
7 4 2 

Let us perform the following transformation of the dependent and independent variables 
then Eq. © turns into 
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Let us consider now the following factorization scheme [3, ^ . The nonlinear second order equation 

y" + .f{y)y'^ + g{y)y' + h{y) = o, 

where 2/' = = D^iy, can be factorized in the form 



under the conditions 



If we assume [D,, 



[Dv - My)y' ~ My)] [A, - My)] y = Q, 
f{y)^-M 

g{y) = 0i(?!)3y -(l)2~h - -p-y^ 

ay 

h{y) = (t)2(f>3y- 

^3(2/)] y — ^{y)j then the factorized Eq. can be rewritten as 



y' - M = ^, 

biy' + cl)2)n = 0. 



(13) 

(14) 

(15) 
(16) 
(17) 



(18) 
(19) 



We can introduce the functions by comparing Eqs. (112]) and p3)) . Then, 0i = 



^ and (/)3 = 

2ai{D + Eyi~'^), where 01(7^ 0) is an arbitrary constant, are proposed. 

Eq. (IT^ can be easily solved for the chosen factorizing functions obtaining as result ft — where /Ci is 

an integration constant. Then, Eq. (jl8p turns into the equation 



y' - 2a, (d + Ey^-') y - n^e'^l'^^y^l'' = 0, 



whose solution is also solution of Eq. (IT^ . 

Furthermore, the following relationship is obtained from Eq. ([16 



AaiD - a^^ - 2aiD + aiE{A - 3 + 2s)y5-^ = 3 + Cy^ 



(20) 



(21) 



Eq. (I^Tj) is a noteworthy result which provides the explicit form of a, and the relationship among the parameters 
entering Eq. (I12[) . Then, the viscous parameter s as a function of parameter 7 is obtained. By comparing both sides 
of Eq. ((2T|) and assuming r = 1 — -, leads to obtain: 



s(7)d 



±V2- 



v3/2 



273/2 



Then, s- e [0, .25] V7 G [1.2599, 2], and s+ G (.75, 1.2071068] V7 G [1,2). Also, the explicit form of ai is 

27V2 



a (7)1= 



± 



3(V2T7^/')' 



(22) 



(23) 



Then, ai- S [-1/3, -.29499] V7 e [1.2599,2], and ai+ e [1.60947, 00) V7 e [1,2). 
We find the following significative values 



7 


s_ 


ai- 
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The main difference of these results from other approaches is expressed through Eq. ([2^ . which represents an 
advantage of the factorization method as opposed to different approaches studied by other authors. This equation 
provides the relationship between the parameters s and 7 in such a way that by fixing s we get a particular value of 

7- 
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The main dynamical variables of the FE are given in parametric form as follows 

/('?)- /oexp (77- r;o), (24) 

H iv) = y'/^ (V) , (25) 

(26) 

p{v)^3y{v), (27) 

p(r?) = 3(7-l)y(?7), (28) 

Ti{v)^~{siy{v) + ^y (29) 

U'?) = — , (30) 
P 

S (7?) = ~3A;b / ni7j)f{7j)Hirj)T{7j)-'yirj)-'/'drj. (31) 



The authors have not been able to find the most general solution of Eq. (l20l) . However, this equation can be studied 
for some specific cases providing particular solutions of physical interest. In Sections III and IV, the cosmological 
solutions as obtained for the viscosity parameter s ^ 1/2 and s = 1/2 arc studied. 



III. SOLUTION WITH s / 1/2. 



In this section, some particular cases of Eq. ((20)) for s 7^ 1/2 are studied to obtain exact particular solutions of FE 
©-dH). By setting ki = 0, Eq. ^ simplifies as 

y' -2ai(^D + Ey^-')y = 0, (32) 

whose solution is given by 

/ p\2/(2s-l) 
yi^) ^ (^.,6-^(^-1)" - I j , (33) 

where K2 is an integration constant. Therefore, the parametric form of the time function is obtained from Eq. ([TT]) 
as follows 

r r f v\ i/(i-2s) 



A. Case s = 0. 



The first special case considered corresponds to s = 0, which means that the bulk viscosity coefficient ^ — const. 
The following particular solution is obtained 

y{v) = (^^26-'^^^" - , and < (,7) = --L (^26-"^"^ + vEa,) . (35) 

Eqs. and provide the corresponding constant parameters ai_ = —0.295 and 7 = ^^2, respectively. A 

particular equation of state is obtained once again through Eq. ([22]) . In Figs. [T]and[2l the behavior of the FE main 
quantities for different values of constant K2 is plotted. 

As we can see, the solution for K2 = —1 is non-singular since p(0) — const. For K2 — —2 and K2 = —3, the energy 
density has a singular behavior when t = 0, since it runs to infinity when time tends to zero, i.e., p(G) — > 00. The bulk 
viscosity, 11, is negative for all values of t, i.e., 11 (t) < Vt € K"*", which is a thermodynamically consistent result as 
expected for K2 = —1. For K2 = —2 and K2 = —3, the solution is valid only when t > to, i.e., H {t) < \/t > to, while 
n (t — )■ 0) > 0. Then, for this interval of time, t S {0,to), the solution has no physical meaning. The entropy behaves 
like a strictly growing time function; then, there are a large amount of comoving entropy during the expansion of the 
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FIG. 1: Solution with s = 0. Plots of energy density p{t), bulk viscosity Tl{t) and entropy Dashed line for At2 = —1. Solid 

line for K2 — —2. Long dashed line for K2 — —3. 




FIG. 2: Solution with s = 0. Plots of the deceleration parameter q{t) and parameter l{t). Plots of energy density p{t), bulk 
viscosity II(t) and entropy Dashed line for K2 — —1. Solid line for K2 = —2. Long dashed line for K2 = —3. 



universe. The deceleration parameter runs from q{0) = —0.5 to q{t) — —1. Then, the solution is accelerating, i.e., it is 
inflationary. The deceleration parameter tends to — 1 as f — cx) (accelerating solutions) but shows a singular behavior 
when time runs to zero. The parameter l{t) shows that all the plotted solutions are far from equilibrium since they 
are inflationary solutions, which is a consistent result. To the best of our knowledge this solution is new. 



B. Cases = 1/4. 

The second case considered corresponds to s = 1/4. In this case, Eqs. (22) and (|23l) provide ai — —1/3 and 7 = 2. 
Therefore, Eq. (PO)) simplifies as 

y' + 2(3)3/4^5/4 _ ^^e-^^yS/^ = 0. (36) 
If we perform the transformation z = y^/^ in Eq. ()36p . then we get the Riccati equation 



33/4 , 1 



-z^ Kie 



= 0, (37) 



2 4 

whose general solution is given in terms of Bessel J„ and Neumman iV„ functions, 

where C(?7) = f^e'^''/^ and K2 is an integration constant. Therefore, the following special solution for Eq. (p6| is 
obtained: 

y^^^ - [^^^hoiav)) + -2Nom)) ) ■ [^^^^ Mav)) + -2Nom)) ) ^- ^ ^ 

In order to study the behavior of the EE dynamical variables in their parametric form, the calculation of Eq. ([M)) 
has been numerically addressed. The solution depends strongly on the value of the numerical constants, in such a 
way that our solution is physical only for K2 < and for negative and relatively small values (< 20) of ki. Numerical 
analysis of the solution plotted in Eig. [3] shows that the solution is singular since the energy density tends to infinity 
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when t —>■ 0. The bulk viscosity is positive, 11 > 0, in the region (0, <*) so the solution has physical meaning only when 
t > t^:, for this era 11 becomes negative as expected from the thermodynamical point of view and tending to zero in 
the large time limit. In the same interval of time (0,i,) the entropy production is negative, < (unphysical 
situation), nevertheless when t > t^, a large amount of comoving entropy is produced during the expansion of the 
universe. 
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FIG. 3; Solution with s = 1/4 and "f — 2. Plots of energy density p{t), bulk viscosity Il(t) and entropy Dashed line for 

Ki — 4, K2 = —10. Solid line for ki — 19, K2 = —3. Long dashed line for ki = 3, H2 ~ —0.7. 

Regarding the dynamical behavior of solution ([5^ . in Fig. 2] the behavior of parameters q and I has been plotted. 
As we can see, the deceleration parameter shows that the universe starts in a non-inflationary phase, but quickly 
entering a inflationary one since q < 0. The plots of l{t) are consistent with this behavior, showing that the solution 
starts in a thermodynamical equilibrium but in a finite time they are far from equilibrium since they are inflationary 
solutions. 
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FIG. 4: Solution with s = 1/4 and j — 2. Plots of the deceleration parameter q{t) and parameter l{t). Dashed line for m = 4, 
K2 = —10. Solid line for ki — 19, K2 = —3. Long dashed line for ki — 3, K2 — —0.7. 

A similar solution has been obtained by Mak et al [1^ but, as we have shown, our solution is qualitatively different, 
with a very different physical meaning. 



1. A particular solution for the case s — 1/4. 
If we set Ki = in Eq. (l36l) . then we get the very simple ODE 

whose solution is given as 

/ (3)3/4 



yiv) 



and t (77) — -Viif + -3''?y^K2 + Vi'^l^ 



(40) 



(41) 



where K2 is an integration constant. In Figs. [S] and [S] the behavior of the FE main quantities has been plotted for 
different values of the constant K2 ■ 
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FIG. 5: Particular solution for s = 1/4. Plots of energy density p(t), bulk viscosity Tl{t) and entropy Dashed line for 

K2 = 0. Long dashed line for K2 — I. Solid line for K2 = 2. 
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FIG. 6; Particular solution for s = 1/4. Plots of the deceleration parameter q{t) and parameter l{t). Dashed line for k,2 — 0. 
Long dashed line for K2 ~ 1. Solid line for K2 = 2. 



The solution has been plotted for three different values of constant K2- The energy density presents a singular 
behavior only for K2 = 0, while the other two solutions show a non-singular behavior when t = 0. The solution for 
K2 = 2 runs quickly to zero. The bulk viscosity is always a negative time function for K2 = 1 and K2 — 2, but the 
solution for K2 = is valid only for t > to since n(t — > 0) > 0, which means that it lacks of physical meaning in 
the interval of time t S {0,to). The entropy always behaves like a growing time function but for the case K2 — the 
universe starts with a non- vanishing entropy, i.e., S(0) = const., while for the other two solutions S(0) — > 0. The 
plots in Fig. [5] show that a large amount of entropy is produced during the expansion of the universe. Regarding the 
deceleration parameter, the plotted solutions run to an acceleration region since q(t) — —1/2 in a finite time. For this 
reason, the solution starts in an equilibrium regimen but quickly run to a non-equilibrium state as shown by plots of 
l{t). A particular solution of this case has been studied by Harko et al [13] obtaining different behavior of the FE 
main quantities. 



Case s 



The second important case considered corresponds to s = 1. According to Eqs. (|22|) and ([23)l . this solution is valid 
only for the equation of state with 7 = 1.25992. Other authors have already studied similar cases for s = 1, 

but with different equation of state (see for instance [l^l with 7 = 2) obtaining different results. Then, according to 
Eqs. ([55]) and ([M)) . the following particular parametric solution is obtained: 



yiv) 



K2e 



aiDt] 



E 
D 



and 



tiv) 



1 



In 



E 

1^2 



rjDai 



(42) 



Then, the FE main dynamical variables can be explicitly obtained through Eqs. (I24|) -(|3T |) . 

In Figs. [7]and|Sl the behavior of the main quantities by giving different values to the constant K2 has been plotted. 

As we can see the solution is valid only for t > to. The energy density is a decreasing function, but the function 
behaves like a constant for a t > tc. The behavior of the bulk viscous parameter shows that the solution is valid only 
for t > to since the solution is positive when t — >■ 0, decreasing and going to a negative constant value during the 
cosmological evolution, which is consistent from the thermodynamical point of view. In the same way, the entropy 
behaves like a growing function only for t > to, showing that a large amount of comoving entropy is produced. 
Nevertheless, the deceleration parameter shows that the universe starts in a non-inflationary phase, but quickly 
entering a inflationary one since q — >■ — 1 Vk2. The plots of l{t) show that plotted solutions are far from equilibrium 
since they are inflationary solutions. 
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FIG. 7: Solution with s = 1. Plots of energy density p{t), bulk viscosity Il{t) and entropy Dashed line for K2 = 0.1. Long 

dashed line for K2 = 10. Solid line for K2 = 100. 
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FIG. 8: Solution with s = 1. Plots of the deceleration parameter q{t) and parameter l{t). Dashed line for K2 — 0.1. Long dashed 
line for K2 — 10. Solid line for K2 = 100. 



IV. SOLUTION WITH s = 1/2. 



We consider now the very special case s = 1/2. This has been the most important and studied case (see for example 
example [3 , [I^ , [1^ , [HI ) within the framework of the bulk viscous cosmological models since, as it has been pointed 
out for several authors, this solution is stable from the dynamical systems point of view [26| as well as from the 
renormalization group approach (27| . 

In this case, Eq. © reduces to: 



H 

H - Ai— + {3 + Ci)HH + {Di + El) = 0, 
n 



(43) 



where Ai = (1 + r) = 2 - i, Ci = \/3, Di = | (7 - 2), Ei = |\/37, and r = 1 - I/7. Since the coordinate 
transformation given by Eq. (jlip leads to obtain several unphysical solutions for s = 1/2, we perform the more 
suitable change of variables given as follows (see also [18]), 



drj = 3[1 + -j=] Hdt. 



(44) 



Then, Eq. (gSj) turns into 



where — (7 -l- 6) — 
s = 1/2. 



Eq. 



y"-^y'^ + y' + 2jby = 0, (45) 
5|) can be solved by factorization providing new exact parametric solutions for 



Eq. (|45p admits the factorization 



2y ' , 



[D - 2ai76] 2/ = 0, 



(46) 
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which can be rewritten in the form 

V' - 2a,jbij = (47) 
" {i/y " " " °' ^^^^ 

or equivalently, 

y' -2aijby-kie'^/'''y^^'^ ^0, (49) 
where ki is an integration constant, with solution given as 

'<"^-*(^^^^-')^ 

where Ci is an integration constant, and the parameter oi is restricted to values given by 



4^37 ± J72 (72^/3 - 6O) -973 + 7 (432n/3 - 756) 

«i± = 7 7= ^ 1 (51) 

3(7-4^3 + 6) 

i.e. , ai+ e [—64.31,-8.38], and ai_ G [-0.23,-0.01]. In the following Subsections IV. A and IV. B several possible 
cases of interest are studied. 

A. General solution. 

In this case it is possible to find a explicit parametric equation for t (from Eq. (l50|) ) with Ci ^ 0. It is given as 
follows 



3 



aiki / / —27^ 



^ 272 giexp(^(27-«ij?)) \ 
' aiB — 27 ' aiki I 



To the best of our knowledge the solution given by Eqs. ([50]) and ([5^ has not been previously reported. 
The FE main dynamical variables are given in parametric form as follows 



(52) 



/(»?) = /oexp(77-77o), (53) 
i/ (77)= 2/1/2 (^)^ (54) 

(2 + B) Ci + 2 (ai + 7) ki exp (27 - aiB)) 

liv) = 7 7 TT ^, (55) 

2 [Ci + aiki exp (27 - aiB)J J 

P (»?) = 3?/ (77) , (56) 
p(r;) = 3(7-l)2/(77), (57) 

Ci(i? + 37) + (2 + 3ai) 7ki exp (27 - aiB)) 

n iv) = J W2_ Ly (^) (58) 

Ci+aikiexp(2^(27-aiB)j 

i:{r,)^^e^^{3yf\ (59) 

n(?7) 



(60) 



In Figs. |9]and[T0l the behavior of the FE main quantities has been plotted. The following constant values have been 
chosen: ai+ as given in Eq. dCT]) while B = (7 + 6) - 3ai, ki = 2, Ci = -1, and 7=1, 4/3, 2 as usual. The 
solutions with ai_ are unphysical. 
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FIG. 9: Solution for s — 1/2. Plots of energy density p{t), bulk viscosity Il{t) and entropy T,{t). Solid line for 7 = 2. Long 
dashed line for 7 — 4/3. Dashed line for 7=1. 
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FIG. 10: Solution for s = 1/2. Plots of the deceleration parameter q{t) and parameter l{t). Solid line for 7 = 2. Long dashed 
line for 7 — 4/3. Dashed line for 7 = 1. 



The energy density shows a singular behavior as t — > 0, but in a finite time it behaves as a decreasing time function. 
This solution is valid for all values of time and 7. The bulk viscous pressure, IT, is a negative decreasing time function 
during the cosmological evolution, IT < Vt G M^, as it is expected from a thermodynamical point of view. The 
viscous pressure also evolves from a singular era but it quickly tends to zero, i.e., in the large limit the viscous pressure 
vanishes as the viscous coefficient, which also becomes negligible small. The comoving entropy behaves as a growing 
time function. There exists a fast growth of entropy for 7 = 4/3, while for 7 = 1 the entropy grows slowly. The 
entropy evolves from a non-singular state, i.e., S(0) — 0, but it quickly grows in such a way that a large amount 
of entropy is produced during the cosmological evolution. The picture of parameter q{t) shows that all the plotted 
solutions start in a non- inflationary phase, but they quickly run to an inflationary era since this quantity runs to — 1 
for all the equations of state. For this reason, the parameter l{t) shows that the solutions are far from equilibrium 
since they are inflationary solutions. 



B. Particular solution 



In the case, it is possible to find a particular solution for t from Eq. ([50]) with Ci — 0. For this case, the solution 
simplifies as follows 



y (77) = cxp {Bi]) 



'aikiexp (2^(27-015)^ 
2-f-aiB 



27 



and t{fj) = (^V3-3)^y-'/^. (61) 
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Then, the FE main quantities are given in the following form: 



fiv) = /oexp (77-770) , 

I N (ai +7) 

P{^) = 3y (ry) , 
p (??)== 3 (7 - 1) y {if) , 
(2 + 3ai)7 



(2 + 3ai)7 



y(?7)' 



(62) 
(63) 

(64) 

(65) 
(66) 

(67) 

(68) 
(69) 



n(77) 

^^^^^-3 a, (7-1) 
E(,7)=7e^"(32/(r/))^/^ 

It is possible to recover the known scaling solution studied by several authors [Hl,!!^ and [sB] from Eqs. (f62|) - ([69l) : 

/ = /oi''° (70) 

H{t)^Hot-\ (71) 

q{t) = H-^ - 1, (72) 

p(^)-Po^-^ (73) 

p (t) = 3 (7 - 1) pot-\lV (t) = -Hop (t) , (74) 
Ho 



lit) 
E(i) 



3(7-1)' 

7^0 
37iIo - 2 



(75) 

(76) 



where Ho = , , fc„^ = 1, Eq 



-3n, 



> 0, and Ho > 0. 



— 67 1,-1 

"_(N/3-3)ai ' 

In Figs. [TTjand[T2l the FE main quantities have been plotted using the same numerical values of Figs. l9land[T0l 
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FIG. 11: Solution for s = 1/2 and Ci = 0. Plots of energy density p(t), bulk viscosity \l{t) and entropy S(t). Solid line for 
7 = 2. Long dashed line for 7 = 4/3. Dashed line for 7 = 1. 

All the plotted solutions have physical meaning V t. These solutions and the ones presented in the last solution (with 
Ci ^ 0) have a similar behavior. We get the following numerical values for parameters q{t) and l{t): qi — —0.98569, 
54/3 = -0.92081, q2 = -0.796 58, while I4/3 = 3.84 and h = 1.86. 



V. SOLUTIONS THROUGH THE LIE GROUP METHOD FOR s = 1/2 



In order to find new solutions and compare with the ones obtained through the factorization method, we study the 
Hubble rate Eq. (gS]) with s = 1/2 through Eq. (gS]) by applying the Lie group method [3l|. Eq. (gSJ admits the 
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FIG. 12: Solution for s = 1/2 and Ci = 0. Plots of the deceleration parameter q{t) and parameter l{t). Solid line for 7 = 2. 
Long dashed line for 7 — 4/3. Dashed line for 7 = 1. 



following symmetries: 



a = [1,0], 6 = [0,y], 6 = [1,2/], 

^4,5 = 0,2/^exp(|(±a-l)^ 



^6 



(1-4) 



exp (I (1 T a; 



1 ± a 

A- 2 



(1-4) 



exp (I (1 T a] 



(77) 



where a = \/T~~8B~\~4AB , and S = (7 + 6) — |. The non-zero constants C^,j verifying the relationship [a,^] 



C^i^k are 



[6 , '^4] = (^14 C4 , [Cl , '^5] — Cfs^ , [^ , a] - €'24^,4 , [6 , Cs] — (^25^5 



(78) 



Then, we shall try to find a suitable change of variables with the symmetries ^ and S,^. These symmetries, = 
[0, y'*/2e''/2(='=''~i)] , bring us to get the following cv that will transform the original ODE into a quadrature. Following 
the the standard procedure we get: 



which lead us to obtain the following ODE and the corresponding solution: 

Then, the solution to Eq. P5)) is given as follows 

V 27a J 



27 



(79) 



(80) 



(81) 



where a = y/1 + ABA - 8B, A = 2- - and i? = ^ (7 + 6) - |. In the following Subsections V.A-V.D, the solutions 



provided in Eq. ([81]) are separately studied. 



For the solution 



A. Solution y- with C2 / 



y- in) = f--^Cie^''(-«-i) + C2e^''('^-i: 
V 27a 



27 



(82) 



with C2 7^ 0, it is possible to find an explicit parametric equation for t through Eq. (|44|) . It is given as 



iv) = 



(3 - V3) (1 - ^^n22^ 



3(1 + a) 7?/ 



TJT 



2F1 



7, ■ 



(l + a)7 7 + a(2 + 7) 2a-/C2e'"i 



2a 



2a 



(83) 
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As we can see, a similar solution to the one obtained through the factorization method has been found. However, as 
it is shown below, they present several important differences. 

The FE main dynamical variables are given in parametric form as follows 



fiv) = /oexp (?7-?7o), 

Hiv) = y 



piv) 
p{v) 

n(r;) 
E(ry) 

liv)) 



2Q7C2e°'' (2 + 7 (q - 1)) + Ci (7 (a + 1) - 2) 
2a7C2e"'') ' 



2(Ci 

Sy- (??) , 

3(7- l)y- {i]), 

7 (2a7 (g + 2) Cae"' + Ci (a - 2)) 

Ci - 2a7C2e'''' 
76^" (3z/_ (r;))!/-^ , 
|n(r7)| 



y- iv) , 



p{i]) 

In Figs. [T51 and [TH the behavior of the FE main quantities has been plotted. 
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FIG. 13: Solution for i/_ with C2 7^ 0. Plots of energy density p{t), bulk viscosity Il{t) and entropy T,{t). Solid line for 'y — 2. 
Long dashed line for 7 — 4/3. Dashed line for 7=1. 
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FIG. 14: Solution for y_ with C2 7^ 0. Plots of the deceleration parameter q{t) and parameter l{t). Solid line for 7 = 2. Long 
dashed line for 7 = 4/3. Dashed line for 7 = 1. 



As it is shown in Fig. 1131 the solution is not valid for 7 = 4/3. For 7=1 (matter predominance) and 7 = 2 
(ultra-stiff matter), the energy density behaves as a decreasing time function during the cosmological evolution. This 
solution is valid for all values of time, except in the case 7 = 4/3, where < 0. The bulk viscosity is a negative 
increasing time function, except in the case 7 = 4/3, where 114/3 > 0- The energy-density, bulk viscosity and entropy 
have a very similar behavior for the cases 7 = 1 and 7 = 2. The solution has a singular origin since the energy density 
tends to infinity as i — >■ 0. The entropy is a growing time function which shows a large amount of comoving entropy 
during the expansion of the universe. In the case 7 = 4/3, the entropy starts growing at t = 60, although we have ruled 
out this case. The behavior of parameter q(t) shows that the solution for 7 = 2 starts in a non-inflationary phase, 
but after a period of time the solution enters an inflationary era. Nevertheless, the solution for 7 = 1 is inflationary 
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for all values of t. The behavior of parameter l{t) shows that the solution for 7 = 2 is close to equilibrium, which is 
thermodynamically consistent. 



B. Solution y_ with C2 = 



For the case j/_ with C2 ~ 0, the solution is given by (after simplifying) 

27 



and 



t- iv) 



(2/- 



-1/2 



drj 



(3-%/3) _i 



3(1 + a) 7 



y- 



/2 



and the FE main dynamical variables are given in parametric form as follows 

fiv) = /oexp(r/-ryo), 
H (77) = yi^^ (ry) , 

q{v) = i((l + a)7-2), 

P iv) = 3y_ (r/) , 

p (r?) = 3 (7 - 1) y_ (7?) , 

Ui-n) = (a-2)72;_ (r?) , 



(«-2)7 



E (7?) = 7e'" (3y- (77))^/^ . 



(92) 

(93) 
(94) 

(95) 

(96) 
(97) 
(98) 

(99) 
(100) 



The behavior of the FE main quantities has been plotted in Figs. [TSjandfTH] As it is observed, in this case, we may 
recover the scaling solution. 
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FIG. 15: Solution j/_ with C2 = 0. Plots of energy density p{t), bulk viscosity Il{t) and entropy S(f). Solid line for 'y — 2. 
Long dashed line for 7 = 4/3. Dashed line for 7=1. 



In this case, as in the last solution with C2 7^ 0, the solution for 7 — 4/3 is unphysical. All the main quantities 
behave in the same way as the last solution with C2 ^ described above. Nevertheless, it is found that qi = —0.015 
for 7 = 1 which represents an inflationary solution, and (72 = 0.732 for 7 = 2 which represents a non-inflationary 
behavior, while I2 — 0.845 < 1, i.e., the solution is within an equilibrium regime. As it has been shown, most of 
the known exact solutions of the gravitational FE with a viscous fluid do not satisfy the condition Z < 1, i.e., the 
condition of thermodynamic consistency, since they show an inflationary behavior. In the case for 7 = 2, we have 
obtained a solution which is thermodynamically consistent and it may describe the early dynamics of a super-dense 
post-inflationary era when the dissipative effects produced by the bulk viscosity may play an important role. 



C. Solution y+ with C2 7^ 



For the solution given by 

y+ iv) = (^C'ie^''('^-i) + C2e-^''('^+i)^ ^ , (101) 
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FIG. 16: Solution y- with C2 = 0. Plots of the deceleration parameter q{t) and parameter l(t). Solid line for 7 = 2. Long 
dashed line for 7 = 4/3. Dashed line for 7 = 1. 



with C2 7^ 0, we get the explicit parametric equation for the time function 



t+ ('?) 



(3-V3) + 



7, 



(l + a)7 7 + a(2 + 7) Cie'"' 



2a 



2a 



3 (l + a)7//2 

The main dynamical variables of the FE are given in parametric form as follows 

f+{ri) = /oexp {v~Vo), 

H+ iv) = iv) , 

^ 2a-fC2 (7 (a + 1) - 2) - Cie"' (7 (a - 1) + 2) 
2(^1-207^26'^'') 

P+ iv) = 3y+ iv) , 

P+ iv) = 3 (7 - 1) y+ (77) , 

7 (2a7 (a - 2) C2 - Ci (a + 2) e"*") 



207^^2 



n+ (v) 



y+, 



Cie^"' + 2076*2 
7 (2a7 (a - 2) C2 - Ci (a + 2) e"') 



(7 - 1) (Cie'") + 207^2) 
E+ (,7) = 76^" (3?/+)'/" 
We have plotted the behavior of the FE main quantities in Figs. [TTlandlTSl 
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FIG. 17: Solution y+ with C2 7^ 0. Plots of energy density p{t), bulk viscosity Il{t) and entropy S(f). Solid line for 7 = 2. 
Long dashed line for 7 = 4/3. Dashed line for 7 = 1. 



This solution shows a behavior quite similar to the one obtained through the factorization method. The energy 
density is a decreasing time function and it is valid for all values of time. The bulk viscous pressure is a negative 
increasing time function, while the entropy is a positive growing time function. As in the case of the factorization 
method, the obtained solution is valid for all the possible values of parameter 7. We find a fast growth of entropy for 
7 = 2, while it grows slowly during the evolution of the universe for 7 = 1. The behavior of parameter q{t) shows 
that all the plotted solutions start in an inflationary phase, since this quantity is close to —1 for every value of 7. The 
behavior of parameter l{t) shows that the solutions are far from equilibrium since these are inflationary solutions. 
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FIG. 18: Solution y+ with C'2 7^ 0. Plots of the deceleration parameter q{t) and parameter l{t). Solid line for 7 = 2. Long 
dashed line for 7 — 4/3. Dashed line for 7=1. 



D. Solution J/+ with C2 = 

In the case of solution y+ with C2 = we get 

.r^r^,.H-)V^ and .-_(v^-3) 



y+ iv) - l^^^ie^'^""'^ 



The FE main dynamical variables are given in parametric form as follows: 

fiv) = /ocxp(77-77o), 

qiv) 



in) 

1 



P iv) = 3y+ (77) , 

P (?7) = 3 (7 - 1) y+ (77) , 

n(r/) = (a + 2)7y+, 



(a + 2)7 



(7-1) 
S(,7)=7e^'' (3y+)'/''. 



We may recover the scaling solution as above. 
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FIG. 19: Solution y+ with C2 = 0. Plots of energy density p{t), bulk viscosity Ii{t) and entropy S(f). Solid line for 7 = 2. 
Long dashed line for 7 = 4/3. Dashed line for 7 = 1. 



In Figs. [19] and [20l the behavior of the FE main quantities has been plotted. As it can be seen, a very similar 
behavior to the scaling solution obtained through the factorization method has been obtained. Therefore, we get the 
same description and conclusions. It is worth mentioning that the following values for the deceleration parameter 
q{t) are obtained: qi = —0.984206, (74/3 = —0.905604, and q2 = —0.732051, while for parameter l{t) we obtain 
h/3 = 3.81121, and I2 = 1.82137, i.e., the same values as the ones obtained for the scaling solution. 
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FIG. 20: Solution j/+ with C2 ~ 0. Plots of the deceleration parameter q{t) and parameter l{t). Solid line for 7 = 2. Long 
dashed line for 7 — 4/3. Dashed line for 7=1. 



VI. CONCLUSIONS. 



In this work, we have studied a flat FRW cosmological model with a matter model described as a full causal bulk 
viscous fluid. By assuming the state equations given in Eq. ([7]), the cosmological model simplifies to a nonlinear 
second order ODE, the Hubble rate equation, for which a coordinate transformation is performed in order to apply 
the factorization method. Due to the coordinate transformation developed on the Hubble rate equation, parametric 
exact solutions have been found. The standard procedure of factorization provides the first order ODE (PH]) . and the 
restriction condition given in Eq. (I2ip which provides a relationship between the viscous parameter s and 7. Then, 
the analysis developed through factorization allows to study the model for all the values of s determined by Eq. (|22p . 
instead of constructing a particular ODE for a single given value of s and arbitrary or specific values of 7, as it has 
been previously studied by several authors. 

We have studied several models for different values of s. Firstly, we have studied and discussed the model for s = 0, 
and 7 = v^. The second model is studied for s = 1/4, and 7 = 2, finding two solutions. The third case corresponds 
to s = 1, and 7 = \/2 « 1.25992. For the very special case s = 1/2, the restriction equation (|21l) provides the explicit 
form of parameter ai. However, the obtained solutions have not restriction on the values of 7. For this important 
case, we have been able to obtain a new solution which reduces, as particular solution, to the known scaling solution. 
To the best of our knowledge, the parametric solutions obtained for all these cases are new. 

In order to obtain more new solutions, the case s = 1/2 has been studied through the Lie group method. The 
analysis carried out allows to obtain two solutions. The solution (|82]) - ([83|) is new, and solution (|101l) - (|102p presents 
the same behavior as the one obtained through the factorization method. Regarding the solution (|82p - ([55|) . it is 
pointed out that it is not valid for all state equation 7. It has been shown that for 7 = 4/3 the solution is unphysical, 
while for 7 = 2 it is thermodynamically consistent and could be relevant from the cosmological point of view. 
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